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USERS' MANUAL FOR COMPUTER PROGRAM FOR ONE-DIMENSIONAL 


ANALYSIS OF COUPLED-CAVITY TRAVELING WAVE TUBES 
by Thomas A. O'Malley and Denis J. Connolly 
Lewis Research Center 

SUMMARY 

The use of the coupled- cavity traveling wave tube (TWT) for space communications 
has led to an increased interest in the efficiency of these devices. Efficiency improve- 
ments have generally been along two lines; (1) the development of efficient multistage 
depressed collectors to recover a hi^ percentage of the kinetic power in the spent 
beam; and (2) improvements in the efficiency of the basic interaction process through 
velocity resynchronization and other methods. To pursue the second line, we have re- 
cently developed a flexible, large- signal computer program written in FORTRAN IV for 
use on the IBM 360/67 time- sharing system. The present report is a users' manual for 
this program. 

The basic programming approach consists of dividing the beam into disks and cal- 
culating the trajectories of these disks as they pass throu^ a sequence of cavities each 
of which may have different geometrical and electrical properties. Extreme flexibility 
is provided in the variety of tube features that can be modeled since each cavity has in- 
dividually entered input parameters. The program can handle lumped or distributed 
severs, input and output couplers, cavity match details provided for cavities near the 
end of a stack, voltage-jump velocity resynchronization with an arbitrary number of 
discrete steps (up to one per cavity), and velocity taper designs of almost arbitrary 
complexity. Backward waves can be handled by an iterative procedure. 


INTRODUCTION 

The use of the coupled- cavity traveling wave tube (TWT) for space communications 
has led to an increased interest in the efficiency of these devices. Efficiency improve- 
ments have generally been along two lines: (1) the development of efficient multistage 



depressed collecters to recover a hi^ percentage of the kinetic power in the spent 
beam (ref. 1); and (2) improvements in the efficiency of the basic interaction process 
through velocity re synchronization and other methods (ref. 2). The NASA Lewis Re- 
search Center has been interested in the first approach to efficiency enhancement for 
some time (refs. 1 and 3). In order to pursue the second, we have recently developed 
a flexible, large- signal computer program written in FORTRAN IV for use on the IBM 
360/67 time- sharing system (ref, 4). The present report is a users’ manual for this 
program. 

The basic programming approach generally follows that described by Vaughan 
(ref. 5), wherein a disk- model beam is followed through a sequence of cavities each of 
which may have different geometrical and electrical properties. Extreme flexibility is 
provided in the variety of tube features that can be modeled since each cavity has indi- 
vidually entered input parameters. The program can handle lumped or distributed 
severs, input and output couplers, cavity match details provided for cavities near the 
end of a stack, voltage-jump velocity resynchronization with an arbitrary number of 
discrete steps (up to one per cavity), and velocity taper designs of almost arbitrary 
complexity. Backward waves can be handled by an iterative procedure similar to that 
described by Kino, et al. (ref. 6). The program has been used to analyze the beam- 
circuit interaction in the 200-watt TWT used on the Communications Technology Satel- 
lite. The computer results showed excellent agreement with experimental data. 

In this report, we first discuss the computer model of the coupled- cavity TWT. A 
detailed discussion of the program structure follows. The input data, output data, pro- 
gram variables, and program subroutines are then described. Symbols are defined in 
appendix A, and a sample problem is given in appendix B. 


MODEL OF COUPLED- CAVITY TRAVELING WAVE TUBE 

The model of a coupled- cavity TWT is illustrated in figure 1. The beam is treated 
as a series of rigid but mutually penetrable disks. A complete description of the beam 
trajectory is obtained by following the disks contained in a single beam wavelength 
(ref. 7). The number of disks per beam wavelength is an adjustable input parameter 
but is usually set equal to 24, The diameter of a beam disk (2b in fig. 1) is allowed to 
take on as many as three different values during ttie course of a calculation. This op- 
tion is provided to accommodate the fact that TWT beams tend to expand in the satura- 
tion region of the tube. The beam may expand until some portion of it is intercepted by 
the tunnel wall. To accommodate this interception, the program allows the disks to 
lose some portion of their charge and mass at specified points along, the tube. 

The tube body is treated as a conducting tunnel of radius a divided axially into a 
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th 

series of discrete cells, where the length of the cell is denoted by Li,. In the cen- 

th 

ter of each cell is a gap of length Impressed across the k gap is a complex 

voltage 

In the absence of a beam, there will be a simple known relation among the various 
Vj^. For a forward- propagating wave in the pass band of a uniform structure, 


'k+i 




= e 


V,, 


( 1 ) 


where o; and are known from cold-test measurements. The power flowing along 
the structure will be given by 


P 


k 



( 2 ) 


where Zj^ also is known from cold-test measurements. 

The voltage wave in the presence of the beam is calculated by adding a complex in- 
duced voltage to the propagating wave, that is. 
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V, + AV, 


k+1 


( 3 ) 


where is an induced voltage to be defined later. 

The body of a coupled- cavity TWT is supposedly a periodic structure. If this were 
the case, all cavities would have the same L, I, o'L, and Z. However, areal 

tube has input and output couplers, severs, and perhaps velocity tapers and cavity 
match details, all of which require individual nonperiodic treatment in the model. 

Thus, the properties of individual cavities are separately specified input parameters in 
the program, and a great variety of tube design variations such as complicated velocity 
tapers or voltage-jump configurations can be easily modeled. Equation (3) is easily 
modified to handle such a quasi- periodic structure. 


= ^k 





where 
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The factor 




ensures that power flow is conserved. 


PROGRAM STRUCTURE 

The program consists of two major elements: the dynamics of the beam disks and 
the beam-wave coupling analysis. These two elements are described in detail in the 
following sections. 


Dynamics of Beam Disks 

The computational problem discussed in this section is the calculation of beam disk 
trajectories in a suitable manner. Ordinarily, the computational procedure is to first 
solve for the disk acceleration. Numerical integration of this equation then yields 

fVi 

v.(t), the velocity of the r” disk at time t. Numerical integration of v.(t) yields 
z^(t), the axial position of the i disk at time . t. 

Instead of using the functions v.(t) and z.(t) to describe the trajectories, it is more 
convenient for our purposes to let z be the independent variable and to calculate the 
functions dt.(z)/dz and t.(z). Here dt.(z)/dz is the reciprocal of the velocity of the 
i*^” disk when it reaches the axial position z, and t.(z) is the time of arrival of the i“ 
disk at the axial position z. In terms of these new functions, the equations of motion 
are 


^ 3/2 



fli 

where F^. is the total axial force on the i ^ disk, c is the velocity of light, m^ is 
the rest mass of a disk, and N^j is the number of disks. 

Instead of using t. and z we use the normalized variables 0^ and ^ defined by 
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and 



where co is the angular frequency and a is the tunnel radius. These normalized var- 
iables are used in the program. At times, however, greater clarity in the explanation 
of a topic can be achieved by using the unnormalized variables t. and z. At such 
times in this repo;rt, the unnormalized variables are used. In terms of the normalized 
variables, the equations of motion are 
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The initial conditions are 
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i coa 
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where is the axial length of a disk and Uq is the initial beam velocity, given by 
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0 


c 



1 





m^c 


In this equation, Vq is the beam accelerating voltage and q^ is the charge of a disk. 
The coordinate system used in the program has its origin at the beginning of the 

fV» 

first cavity. The k^^ cavity is divided into equal parts of length Az^^, given by 
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The n node in the first cavity, denoted by is defined by 




(6) 


The n^^ node in the cavity, denoted by z^, is given by 


^nk 


k-1 

z 

m=l 


" = •••> k = 2, ..., 


( 7 ) 


where N is the number of cavities. A point that lies midway between two nodes z . 

C 1.^ Qj IV 

and z , is denoted by z i . . Such a point is called the n antinode in cavity k 

(fig. 2). The numerical integration of equation (4) is such that the functions dt/dz and 

t^(z) are evaluated at antinodes and the functions d t./dz are evaluated at nodes. 

At some point in cavity k, we have trajectory data at node n - 1 and antinode 

2 2 

n- 1; that is, the functions d t,(z . )/dz , dt.(z i j/dz, and t.(z i ,) are 

1 X y iv 1 n* 2 > ^ ^ n* 2 9'^ 

2 2 

known. The problem now is to obtain d tj(z^)/dz , dT(z^^i ^'i^^n+J^ k^’ 

We first obtain approximations for dt.(z^j^)/dz and with the equations 
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Knowing, at least approximately, dt.(z^)/dz and , we can evaluate 

d^t.(z ,)/dz^ from equation (4) . We now obtain dt.(z ,i i)/dz and t.(z^,i , ) by the 

1 XlK 1 n+’^jK! 1 n“r2>^ 

equations 
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m=l 


(13) 
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where J, is the Bessel function of order 1, is the m zero of the Bessel function 

-L lU 

of order zero, is the permittivity of free space, and and z’ are given by 
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The function F defined in equations (12) and (13) is continuous for all values of z. In 
particular, there is no jump discontinuity at z = 0. 

Using equations (12) and (13) for calculating space charge forces on each integra- 
tion step would be intolerably slow. The method used herein is to set up tables 
corresponding to the discrete disk diameters. Each table contains 100 data points 
corresponding to 100 separation distances. The n*^ data point would be the space 
charge force between two disks, as calculated by equation (12) or (13), separated by the 
n*"^ separation distance. In the computer simulation, the space charge force on one 
disk due to another is calculated by using linear interpolation on the appropriate table. 

If the two disks in question have different diameters, two linear interpolations are done 
on the appropriate tables, and the average is taken. 

Fortunately, the space charge force on a given disk is a slowly varying function of 
the independent variable z. This is especially true in the beginning cavities of the 
tube. To take advantage of the slowly varying nature of the space charge force, the 
calculation is not done on every integration step. The calculation of space charge 
forces is initially done every step, where N is an input. The value of N may 
change during the simulation in a manner to be described later. On those integration 
steps where space charge forces are not calculated, they are. approximated by quadratic 
curve fits. Each disk has its own curve fit. These curve fits are obtained by fitting a 
quadratic polynomial through the latest three calculations of space charge forces. Let 
us say that these three calculations were done when z was equal to z^, z^, ^d z^. 

At an integration step where a new space charge force calculation is due, at z = z^, the 
values calculated are compared with those obtained when the curve fits based on the 
points z^, z^, and z^ are evaluated at z^. If the difference in the values for any disk 
is exceeded by a specified tolerance, the next calculation of space charge forces occurs 
steps later. In general, N is halved each time the specified tolerance is exceeded. 

For calculating space charge forces, the positions of all disks must be known at a 
time t. This information is not readily available since z rather than t is the inde- 
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,th 


pendent variable. To calculate the space charge force on the i “ disk located at some 
z, we must determine the positions of all the other disks at the time t.(z). Consider 
The time of arrival t. (z) at z is known. By the periodicity of the mo- 
at 


the j disk. 

J fVi 

tion, the disk that originally was m wavelengths behind the j disk arrives at z 


the time t.(z) + mT, where T is the reciprocal of frequency. Similarly, the disk that 
J f 

originally was m wavelengths ahead of the j disk arrives at z at the time t (z) + 
(-m)T. We consider all such disks and find m such that 


It^(z) - (tj(z) + mi;) I = Minimum 


(14) 


th til 

We may call the disk in question the j disk in the m*^ cycle. If the input parameter 

th 

KCYCLE is zero, we assume that the j disk in any other cycle is too far away from 
til 

the i disk, at the time t.(z), to contribute significantly to the space charge force on 

th ^ th 

the i*"” disk. If KCYCLE is not zero, we also consider the j disk in cycles m - 1 

and m + 1. Let zP be the position of the disk in the m*"^ cycle at the time t. (z). 

m ^ 

We obtain z, by assuming that this disk has constant velocity v.(z) in the time inter- 


tm 
y t.(z) and t(z) + mT. Then z^ is given by 


J 


m 


= z + Vj(z) [t.(z) - (L(z) + mT)J 


(15) 


Similar expressions are obtained, if required, for The space charge 

force is then calculated based on the separation distance z - z^. The calculation is 

th ^ 

done for all disks, excluding the i disk, and the results are summed to obtain the 

th m 

total space charge force on the i‘'” disk. When the separation distance z - z™ is 

greater than a specified value, the space charge calculation is not done for this term 

in the summation. 

Radiofrequency electric field forces. - Forces on a disk are always evaluated when 
the disk is at a node. This follows from equation (4) and the fact that d t./dz is al- 

fVi ^ fV\ 

ways evaluated at a node. So consider a disk in the k ^ cavity located at the node 
z^. The axial force on the disk due to the radiofrequency (rf) electric field, at time t, 
is given by (ref. 4, eq. (20)) 


^rf= 
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where 
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(18) 


(19) 


(20) 


In equation (16), is the forward- traveling voltage and is the backward- 

traveling voltage for the k^^ cavity. In equations (16) and (17), |ii is a shaping factor 
that is chosen to best suit the geometry of the gap. In the case of very blunt tunnel tips, 

fV» 

11 = 0 . In equation (18), is the lowest order phase shift for the k cavity. In 
equation (20), y^ is the axial position, relative to the center of the k^^ cavity, of the 
n^^ node in the cavity. 

We define by 




"^^mk^nk 


n = 1, . . . , k = 1, . . . , 


,th 


Then the electric field force on the i disk located at is given by 

io;ti(z^)- 


^ink=I^" 


Wlk-^^nk^bk/" 


( 21 ) 


( 22 ) 
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The are calculated and stored in a table when cavity k is entered. The number 

of terms included in the summation is determined by an input parameter. If the number 
of terms in the summation is no larger than about 40, the computer time for calculating 
the table is small compared with the total computer time. If cavity k has the same 
Lj^, /Sqj^, and b^ as cavity k - 1, then and a new table does not have 

to be calculated. When a new table does have to be calculated, the old table is no longer 
needed, and the new table may occupy the same storage space as the old table. The 


electric field force is evaluated by equation (22). 

th 

Voltage- jump electric field forces. - Consider a disk in the k*^” cavity located at 
the n node z^. Let Vj^ be the voltage jump associated with cavity k. Then the 
electric field force on the disk is given by (ref. 4, eq. (31)) 


^nk ^nk^Jk 


where is given by 
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where k is given by 

m 



( 25 ) 


The are calculated and stored in a table and have the same features as the 

for rf electric field forces. 


Beam- Wave Coupling Analysis 

The parameters of interest in the beam- wave interaction in the k*'” cavity are the 
forward induced voltage and the backward induced voltage AV^^j^. First, we de- 

fine the fundamental Fourier component of the current density by 




dt' 


It can be shown (ref. 4, eq. (81)) that Jj^(z) is given by 



i=l 


where is the number of disks. The induced voltages for the k cavity are then 
given by (ref. 4, eqs. (75) to (77)) 
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where is given by 


\ 

In these equations, the integration variable y is the axial position relative to the center 
of cavity k. The integrals in equations (27) and (28) are evaluated using Simpson's rule 
with end corrections. The end corrections were derived with the assumption that the 
number of nodes per cavity N^, is a multiple of four. The results are 


N„ 
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n=l m=-oo 
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where, from Simpson's rule with end corrections, the 6^^ are 


6, =l4-iZ 

^ 3 24 


6 = — , 4 < n < N -2 , if n is odd 
"3 ^ 
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3 24 
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3 24 
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We define by 






‘^mk^nk 


m=- oo 


( 32 ) 


The are calculated and stored in a table that has the same features as the 

and tables. Our final expressions for AV^j^ and are 


N. 




n=l 




^^bk ^ ^n^nk'^l^^nk^ 


n=l 


(33) 


(34) 


Calculating the effects of the backward wave requires an iterative procedure. In 
the first pass throu^ the tube, the AV^j^ are calculated and stored for each cavity k. 
From the AV^j^, the backward voltages for a cavity chain from cavity kl to cavity k2 
are obtained from 


,k2 ^^b.k2 


(35) 




-i(/?iL)k 


k = k2 - 1 , . . . , kl 


(36) 


By a cavity chain, we mean a sequence of consecutive cavities whose backward voltages 
are to be calculated. After are known for the desired cavities, a second pass is 

made through the tube. The second pass yields a new set of AV^j^ that can be used for 
calculating a set of Vbk for a third pass. The process continues until convergence is 
obtained. In many applications the iterative procedure will converge faster if the itera- 
tion is done on one cavity chain at a time. This would be true, for example, for a tube 
having severs. The program has the capability of performing the iterative procedure 
for an arbitrary set of cavity chains. 

The computational procedure in the beam-wave interaction process is as follows: 
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(1) When cavity k is entered, calculate and store the summation in equa- 
tion (32) for each node n. If cavity k has the same L^, and as cavity 

k - 1, then k 1 ^ table does not have to be calculated. 

(2) Obtain a first approximation to AV£j^ by assuming that the disks have constant 
velocity in cavity k. Thus, 




dt.(z^k) 


dz 


^^nk " ^Ik^ 


n = 1, 


N„ 


(37) 


Kno\ving tj^(z^), we can calculate from equation (26) and then calculate AV^j^ 

from equation (33) . 

(3) Let the forward voltage for cavity k be 




-(Q-L) 

e 


f,k-l 


-i(^lL)k_i 




(38) 


If there is a backward voltage, we assume it is known from equation (36) using the set 
of AV^j^ obtained from the previous pass through the tube. 

(4) Proceed with numerical integration of the equations of motion to obtain trajec- 
tory data at each node and antinode in the cavity. 

(5) At the n*"^ node in the cavity, after t.(z , ) has been obtained, calculate and 

* ^ th 

store ^n^nk'^l^^nk^ ^n^nk”^! ^^nk^ ' These quantities are the n terms in the sum- 

mations of equations (33) and (34) . 

(6) When the last integration step in cavity k has been done, calculate a better 
approximation to AV£j^ from equation (33). The AV^^ so obtained replaces the old 
AVfk calculated in step 2 , Also, calculate AV^^^ from equation (34) . The AV^j^ are 
stored for use in the next pass through the tube. 

(7) If additional accuracy is required, make another pass through the cavity, re- 
peating steps 3 to 6. 

(8) Repeat steps 1 to 7 for cavity k + 1. 

(9) When the last cavity is done and the effects of the backward wave have to be de- 
termined, make a second pass through the tube. The backward voltages for the second 
pass are determined from the AVj^j^ of the first pass. 

(10) Make as many passes throu^ the tube as are required to obtain convergence . 
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Summary of Program Steps 

The entire computational procedure is given in flow chart form in figure 3. The 
procedure is summarized by the following steps: 

(1) Read input data. 

(2) Calculate the tables for space charge forces. 

(3) Begin numerical integration of equations of motion with z as the independent 
variable. 

ttl 

(4) When a new cavity is entered, let us say the cavity, print data for cavity 

k - 1 . If the parameters of cavity k are different from those of cavity k - 1, calcu- 
late the tables that are required for evaluation of rf forces, voltage-jump forces, and 
induced voltages. Obtain a first approximation to AV^j^ by assuming that the disks 
have constant velocity throughout cavity k. Attenuate and phase shift and vec- 

torially add the result to AVjj^ to obtain the forward voltage for cavity k. 

(5) When the last integration step in cavity k has been done, calculate a better 
approximation to AV^j^, and use the new approximation in place of the old. Also, cal- 
culate and store 

(6) If additional accuracy is required, make a second pass through cavity k. 

(7) Repeat from step 4 for cavity k +1. 

(8) If the effects of the backward wave are to be determined, make a second pass 
through the tube. 

(9) Make as many passes throu^ the tube as are required to obtain convergence. 

DESCRIPTION OF INPUT DATA 

The following input data are required by the program; 

Description 

Tunnel radius, cm 

Voltage attenuation of forward wave in going from 
cavity k to cavity k + 1, dB per cavity; k = 1, 
LASTCV 

Voltage attenuation of backward wave in going from 
cavity k + 1 to cavity k, dB per cavity; k = 1, 
LASTCV (ALPHLR(K) needs to be loaded only if 
KLOSS= 1.) 


Name Symbol 

ACM a 

ALPHL(K) (“L)fk 

ALPHLR(K) 
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Name 


Symbol 


Description 


B1LDP(K) 

BCM(I) 

COSHML 


DATTEN(I) 


FREQGH 

lOBMA 

DSriTCV 

ISAVE 


7T 

b 


cosh (fx) 


Phase shift of voltage for cavity k; k = 1, LASTCV 

beam radius, cm; 1 = 1, NUMB (See NUMB, NBl, 
and NB2 for related inputs.) 

Shaping parameter for electric field 

E (a, z) at gap edge 

= cosh (fx) 

E (a, z) at gap center 

til 

In interception region, the charge (mass) of a 
disk is the product of DATTEN(I) and the initial 
charge (mass); 1=1, NUMA. The interception 
region is a region throughout which the charge 
(mass) of a disk is constant. DATTEN(l) is unity 
and must be loaded. (See NUMA, NAl, and NA2 
for related inputs.) 

Frequency, GHz 

Beam current at tube entrance, mA 

Number of first cavity considered in the present case 

ISAVE = 11 +12+13, where II, 12, and 13 are defined 
as follows; 

II = 0: No action 

11 = 1: Initial state of the present case will be 
stored in COMMON/STATE/ for use in a future 
case. 

12 = 0; No action 

12 = 2; Initial state of the present case will be 
stored in a specified data set. 

13 = 0; No action 

13 = 4: Final state of the present case will be 
stored in a specified data set. 
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Name 


Symbol 


Description 


ISTATE Used only when INITCV ^ 1: 

ISTATE = 1: Initial state of present case is equal 
to initial state of previous case. 

ISTATE = 2; Initial state of present case is equal 
to final state of previous case. 

ISTATE = 3: Initial state of present case is equal 
to initial state of the last case for which II (see 
ISAVE) was equal to 1. 

ISTATE = 4: Initial state of present case is equal 
to the state stored in a specified data set. 

JSCF Initially, space charge forces are evaluated once 

every JSCF steps. 

K CYCLE K CYCLE = 0; Consider only one rf cycle in calcu- 

lating space charge forces. 

K CYCLE ^ 0: Consider three rf cycles in calcu- 
lating space charge forces. 

KIMP KIMP = 0: Pierce impedance is entered as input, 

KIMP = 1: Total impedance is entered as input. 

(See ZIMP for related input.) 

KLOSS KLOSS = 0: ALPHLR(K) will be set equal to 

ALPHL(K); thus, ALPHLR(K) does not need to be 
loaded. 

KLOSS = 1: ALPHLR(K) has to be loaded. 

KPLOT KPLOT = 0: Plots are desired. 

KPLOT ^ 0: Plots are not desired. 

KPRINT KPRINT = 0: Print input data. 

KPRINT ^ 0: Do not print input data. 

KREL KREL =0: Use relativistic equations of motion. 

KREL ^ 0: Use nonrelativistic equations of motion. 
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Name 


Symbol 


Description 


KSMSIG KSMSIG = 0: Print small- signal parameters, 

KSMSIG 7 ^ 0; Do not print small- signal parameters. 
KSPACE KSPACE = 0; Calculate space charge forces. 


KSPACE ^ 0; Set space charge forces equal to zero. 


KVEL(I) 


Number of cavity for the printout of disk veloci- 

ties; 1=1, NVEL (See NVEL for related input.) 

LASTCV 


Number of last cavity considered in the present case; 
LASTCV < 650 

LCIRCM(K) 

Lk 

Length of k*"^ cavity, cm; k=l, LASTCV 

LGAPCM(K) 


Length of k^^ gap, cm; k = 1, LASTCV 

NAl 


Charge and mass of a disk change from first value to 

fV\ 

second value at end of NAl cavity (NAl needs to 
be loaded only if NUMA >1.) 

NA2 


Charge and mass of a disk change from second value 

th 

to third value at end of NA2''” cavity (NA2 needs to 
be loaded only if NUMA = 3.) 

NBl 


Beam radius changes from first value to second value 

ii. 

at end of NBl*"” cavity (NBl needs to be loaded only 

if NUMB > 1.) 

NB2 


Beam radius changes from second value to third value 
at end of NB2^^ cavity (NB2 needs to be loaded only 
if NUMB = 3.) 

NBWM 


Number of last cavity considered in calculating back- 
ward wave 

NCAVSS 


Number of cavity for which small- signal parameters 
are calculated 

NDISKS 


Number of disks (NDISKS < 48) 
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Name 


Symbol 


Description 


NMAX 

NUMA 

NUMB 

NVEL 

NXGRID N^ 

NXISC 

PINDBM 

SCMAX 

TOLDV 


TOLSC 


For sums that are stored in tables, excluding space 
charge force tables, NMAX is the upper limit 
(NMAX < 40) 

Number of discrete values for the charge and mass of 
a disk (1 < NUMA < 3) 

Number of discrete values for the beam diameter 
(1 < NUMB < 3) 

Number of cavities for which disk velocities are to be 
printed out (See KVEL for related input.) 

Number of nodes per cavity. NXGRID must be di- 
visible by 4 (8 NXGRID < 64) . 

Number of data points in each space charge force 
table (NXISC < 200) 


Input power, dBm 

If the maximum separation between two disks is 
greater than SCMAX disk lengths, the space charge 
force on one disk due to the other is set to zero. 
SCMAX should equal an integer plus one-half. 


Error criterion for determining whether an additional 
pass through a cavity is required. If 


AVp - AVo , , 

f, new f , old 

V, 

f.new 


> TOLDV 


then an additional pass is made. 


Error criterion for determining whether the frequency 
of calculating space charge forces should be 
doubled. If for any disk the difference between the 
calculated space charge force and the extrapolation 
of the quadratic curve fit is greater than the product 
of TOLSC and the maximum rf electric field force 
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Name 


Symbol 


Description 




in the cavity, the frequency of calculating space 
charge forces is doubled. 

VOB 

Vq 

Beam voltage, V 

VJUMP(K) 

VjF 

The dc voltage jump for cavity, V; k = l. 


LASTCV 

ZIMP(K) 


If KIMP = 0, ZIMP(K) is the Pierce impedance in 
ohms for k*"^ cavity; k = 1, LASTCV. If 
ZIMP = 1, ZIMP(K) is the total impedance in ohms 
for k^^ cavity; k = 1, LASTCV. 


DESCRIPTION OF OUTPUT DATA 

The program output consists of both printed output and plots. The plots are pro- 
duced only if the input parameter KPLOT is zero. The printed output includes three 
parts: (1) printing of small- signal parameters (only if the input parameter KSMSIG is 
zero); (2) a cavity-by- cavity printing of selected data; (3) printing of normalized disk 
velocities at selected cavities. The program can be easily modified to calculate and 
print other data if desired. The format of the output can be seen in the sample problem 
of appendix B. 

First, the input data are printed in NAMELIST format. (The printing of the input 
data can be suppressed by letting KPRINT be nonzero.) Then the small- signal param- 
eters are printed if desired. The small- signal parameters are as follows: 


Name 

Symbol 

Description 

UO 


Initial beam velocity, m/ sec 

BEB 


Product of beam propagation constant and beam 
radius 

BIB 


Product of wave propagation constant and beam 
radius 

KP 

% 

Pierce impedance, ohms 

ZC 

Z 

c 

Total impedance, ohms 


21 



Name 

Symbol 

Description 

C 

C 

Pierce's C 

B 

b 

Pierce's b 

D 

d 

Pierce loss parameter d 

DGAIN 


Asymptotic gain per cavity, dB/ cavily 

QQ 

Q 

Pierce's space charge parameter, QC/C 

A1PA2 

Ai +A2 

Launching loss, dB 

The following data are printed for each cavity; 

Name 


Description 

CAV 


Cavity number 

VMAG 


Magnitude of gap voltage, V 

VPHA 


Phase of gap voltage divided by tt 

ISMAG 


Magnitude of normalized induced current 

ISPHA 


Phase of normalized induced current divided by tt 

VPHA-ISPHA 


VPHA minus ISPHA 

LS GAIN 


Large- signal gain, dB 

POUT 


Output power, W 

EPF 


Efficiency 

PKE 


Change in beam kinetic power referred to initial beam 
kinetic power, W 

PFW 


Power in forward wave, W 

PBW 


Power in backward wave, W 

PLC 


Cumulative power loss, W 
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Name 


Description 


PBAL 


Power balance equal to 

^ I PKE - PJUMP + PFW + PBW + PLC 

‘oVo 

where PJUMP is the power due to all voltage jumps 
up to the current cavity 


ITER Number of iterations of the beam disks traversing the 

cavity 

SC Space charge forces are calculated every SC steps. 

Figure 4 is an example of the plotted output. The figure is a form of Applegate 
diagram. We calculate (p^^, the phase of the forward voltage seen by the i ^ disk 
when it is at the middle of the gap in the k*"^ cavity, as 


<Pki = 


Urn 

L fk J 

Re 

L fk ® J 


>1 






The parameter is plotted as a function of k. There are i N^ plots, one 
every other disk. 


for 


DESCRIPTION OF GLOBAL VARIABLES 

Global variables are variables that are used by more than one subroutine and hence 
are put into COMMON blocks. In the following descriptions, we exclude input data var- 
iables since these have already been described. 


Name 

Symbol 

Description 

A 

a 

Tunnel radius, m 

ALPHLB(K) 


Voltage attenuation of backward wave in going from 
cavity k + 1 to cavity k, dB per cavity; k = I, 
LASTCV 
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Name 

Symbol 

Description 

B(D 

b 

beam radius, m; 1 = 1, NUMB 

C 

c 

Speed of light 

Cl (I) 


th 

Constant associated with r beam diameter; 1=1, 
NUMB 

C2(D 


th 

Constant associated with beam diameter; 1=1, 

NUMB 

C3 


Constant associated with calculation of forces due to 
voltage jumps 

CNORM 

c 

cua 

Normalized speed of light 

COSY (I) 


Array of cosine values used in subroutine TABLE; 
I = 1, 2*NXGRID + 1 

CTHETA 

cos 6 

Cosine of normalized time 

CURRNT 


Fundamental Fourier component of current 

DELV(K) 


th 

Forward induced voltage for k^^ cavity, V; k = l, 
LASTCV 

DELVR(K) 


Backward induced voltage for k^^ cavity, V; k=l, 
LASTCV 

DELX(I) 


th 

Integration step size in cavity string; 1=1, 

NSCAV 

DELXSC 


Spacing between points in space charge force tables 

DPR 


Degrees per radian, DPR = 180 /tt 

DVRTEM(K) 


th 

Backward induced voltage of previous case, for k*^^ 
cavity, V; k=l, LASTCV 

DVTEMP(K) 

^Vfk 

th 

Approximate forward induced voltage for k cavity, 


obtained by assuming disks have constant velocity 
throughout cavity, V; k = 1, LASTCV 
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Name 

ECHARG 

EMASS 

EPSO 

FORIER(K) 

FREQ 

lOB 

ISCAV(K) 

ITER 

ITERSV 

KBEGIN 

KCAV 

KOUNTP 

KPDZC(K) 

KSCAV 

KSCF 

KSTOP 

KX 


Symbol Description 

e Charge of an electron (ECHARG is positive.) 

m^ Mass of an electron 

e 

Cq Permittivity of free space 

Fourier coefficient of fundamental harmonic of in- 
duced current in cavity; k = l, LASTCV 

f Frequency, Hz 

Iq Beam current at tube entrance, A 

til 

Number of cavity string in which cavity is 
located; k = 1, LASTCV 

Present number of iteration of beam disks traversing 
cavity 

Total number of iterations of beam disks traversing 
cavity 

Equals zero at beginning of a case and is set equal 
to 1 immediately thereafter 

Number of current cavity 

Counts points to be plotted 

th 

Ratio of Pierce impedance to total impedance for k 
cavity; k = l, LASTCV 

Number of present cavity string (A cavity string is a 
string of consecutive cavities all having the same 
geometrical and electrical properties.) 

Space charge forces are calculated every KSCF steps 

Equals zero until end of last cavity, when it is set 
equal to 1 

Number of current node, counted from first node of 
first cavity 
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Name 

Symbol 

Description 

KXO 


Initial node of present case, counted from first node 
of first cavity 

KXPRI 


Number of current node, counted from first node of 
current cavity 

KXSC 


Number of node, counted from first node of first 


cavity, at which the most recent calculation of 
space charge forces was done 


LAMB 

X 

Beam wavelength, m 

LCIR(K) 

Lk 

Length of cavity, m; k=l, LASTCV 

LDISK 

^d 

Axial length of a disk, m 

LDISKN 

Z/a 

Normalized axial length of a disk 

LSMALL(K) 


One- half of gap length in cavity, m; k = l, 

LASTCV 

MDISKO) 

^d 

tVt 

Mass of a disk in I interception region, kg; 1 = 1, 
NUMA 

MUA(I) 


til 

r” zero of Bessel function Jq; 1 = 1, 40 

NB 


Number of current beam-diameter region (A beam- 
diameter region is a region throughout which the 
beam diameter is constant.) 

NCASE 


Number of current case 

NCAVPS(I) 


Number of cavities in cavity string; 1 = 1, 

NSCAV 

NCUMCS(I) 


th. 

Number of last cavity in r” cavity string; 1 = 1, 
NSCAV 

NQ 


Number of present interception regicai 

NSCAV 


Total number of cavity strings 

POB 


Linear current density, C/ m 
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Name 


Symbol 


Description 


PBW(K) 

PFW(K) 


th 

Power in backward wave in k''” cavity, W; k = 1, 
LASTCV 

til 

Power in forward wave in k^^ cavity, W; k = 1, 
LASTCV 


PI 

7T 

PI = 3.1415927 

PJUMP(K) 


til 

Power due to all voltage jumps up to k cavity, W; 
k = 1, LASTCV 

PKE(K) 


til 

Kinetic power at end of k*"^ cavity minus kinetic 
power at beginning of tube, W; k = 1, LASTCV 

PLC(K) 


Cumulative power loss at end of k^'^ cavity, W; 
k = 1, LASTCV 

POBAL 


Power balance 

Q DISK (I) 


fVi 

Charge of a disk in r" interception region, C; 
I = 1, NUMA 

SBETOAO) 

/3oa 

t\\ 

Value of for cavities in I " cavity string; 1 = 1 

NSCAV 

SCAl(I) 


I^^ entry in table of space charge forces for first 
beam diameter; 1=1, NXISC+i 

SCA2(I) 


til 

entry in table of space charge forces for second 
beam diameter; 1 = 1, NXISC+l 

SCA3(I) 


I^^ entry in table of space charge forces for third 
beam diameter; 1 = 1, NXISC+i 

SIMPS® 

^i 

I^^ coefficient for Simpson's rule integration; 1=1, 
NXGRID 

SINY® 


Array of sine values used in subroutine TABLE; I = 
2*NXGRID+1 

SLCIR® 


Length of cavities in I^^ cavity string, m; 1 = 1, 
NSCAV 
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Name 

Symbol 

Description 

SLSMAL(I) 


th 

One-half of gap length for cavities in cavity 

string, m; 1 = 1, NSCA V 

SR1MB2 


SR1MB2= ^1- (Uq/c)^ 

SSGAIN(K) 


til 

Small- signal gain at cavity, based on small- 

signal calculations using geometrical and electrical 
properties of cavity number NCAVSS; k = 1, 



LASTCV 

STHETA 

sin 9 

Sine of normalized time 

THT(I) 

6. 

1 

d6. 

th 

Value of normalized time of arrival of r" disk at 
current axial position 1 = 1, NDISKS 

THTP(I) 

1 

First derivative with respect to ^ of 0.; 1 = 1, 
NDISKS 


d^0. 

1 

THTPPa) 

Second derivative with respect to ^ of 0^; 1 = 1, 



NDISKS 

TWOPI 

2ir 

TWOPI = 6.2831853 

TXR(1, 1) 


Real part of I^^ entry in table needed to calculate 
induced voltages; 1 = 1, NXGRID 

TXI(1, 1) 


th 

Imaginary part of r‘ entry in table needed to calcu- 
late induced voltages; 1 = 1, NXGRID 

TXRa, 2) 


th 

Real part of entry in table needed to calculate rf 

electric field forces; 1 = 1, NXGRID 

TXI(I, 2) 


Imaginary part of I^^ entry in table needed to calcu- 
late rf electric field forces; 1 = 1, NXGRID 

TXR(I, 3) 


entry in table needed to calculate voltage-jump 
forces; 1=1, NXGRID 

UOB 


Initial beam velocity, m/sec 
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Name 

Symbol 

Description 

UONORM 

oja 

Normalized initial beam velocity 

V0BSUM(K) 


til 

Beam voltage plus sum of all voltage jumps up to k 
cavity, V; k = 1, LASTCV 

VGAP(K) 


fVi 

Sum of forward voltage and backward voltage for k ^ 
cavity, V; k = 1, LASTCV 

VGAPF(K) 

% 

Forward voltage for k^^ cavity, V; k = 1, LASTCV 

VGAPR(K) 

Vbk 

fVv 

Backward voltage for k cavity, V; k = 1, LASTCV 

VGAPIM 


Voltage across first gap, V 

W 

a; 

Radian frequency, rad/ sec 

X 


Current value of normalized position along tube axis, 
measured from beginning of first cavity 

xo 

^0 

Initial value of normalized position along tube axis 
for the present case, measured from beginning of 
first cavity 

XICAV(K) 


Normalized axial position at end of k^^ cavity; 
k = 1, LASTCV 

XIMXSC 


Maximum normalized separation distance in space 
charge force tables (If two disks are separated by a 
distance greater than XIMXSC, the space charge 
force on one disk due to the other is set to zero.) 

XISCl 


One of three latest values of | at which space charge 
forces were calculated 

XISC2 


One of three latest values of | at which space charge 
forces were calculated 

XISC3 


One of three latest values of ^ at which space charge 


forces were calculated 
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Name 


Symbol 


Description 


XISCAVO) 

XPRI 

XSCA(I) 

XSCBa) 


Normalized axial position at end of r” cavity string; 
1=1, NSCAV 

Present value of normalized position along tube axis, 
measured from beginning of current cavity 

Coefficient in quadratic curve fit of space charge 
force on disk; 1 = 1, NDISKS 

Coefficient in quadratic curve fit of space charge 
force on disk; 1 = 1, NDISKS 


XSCC(I) 

XSCl(I) 

XSC2(1) 

XSC3(I) 

ZC(K) 


Coefficient in quadratic curve fit of space charge 
force on disk; 1=1, NDISKS 

Third most recent evaluation of space charge force on 
disk; 1=1, NDISKS 

Second most recent evaluation of space charge force 
on I^^^ disk; 1 = 1, NDISKS 

til 

Most recent evaluation of space charge force on r” 
disk; I = 1, NDISKS 

Total impedance for k^^ cavity, ohms, k = 1, 
LASTCV 


DESCRIPTION OF SUBROUTINES 

In this section, the COMMON blocks, the main program, and each subroutine are 
described. Flow charts are included for the more complex subroutines. 


COMMON Blocks 

The program has five COMMON blocks. COMMON/INDATA/ contains the input 
data. COMMON/MREAL/ contains real and integer global variables. COMMON/ 
MCOMP/ contains complex global variables. COMMON/PLOT/ contains arrays used 
for plotting Applegate diagrams. COMMON/STATE/ contains arrays used for storing 
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initial states and final states. An initial state is a collection of variables whose values 
are those at the beginning of the first cavity in the tube section under consideration. A 
final state is the same collection of variables but whose values are those at the end of 
the last cavity in the tube section. The variables in this collection are those whose 
values must be known in order to start the simulation of a tube section. It is necessary 
to store these initial and final states because they may be used in later cases. For ex- 
ample, in simulating the second tube section, we use the final state of the first tube 
section as a startup. 


Main Program 

The main program reads the first 40 zeros of the Bessel function Jq. Subroutines 
STDAT, INDAT, and SMSIG are then called. When the last integration step has been 
done, the final state is stored in COMMON/STATE/. Depending on the value of 
ISTATE, the final state may also be stored in a data set for possible use in a future 
case. The main program calls XINT to integrate one step and calls BTWNS to check on 
what action should be taken between steps. 


Subroutine ACC 


Subroutine ACC is called by subroutine XINT. Subroutine ACC calculates the nor- 
malized acceleration, ACCX, of each disk. The accelerations result from three 
forces; space charge force, rf electric field force, and voltage-jump force. The ac- 
celerations due to space charge forces are obtained by quadratic curve fits, and the ac- 
celerations due to rf electric field force and voltage- jump force are calculated in sub- 
routine RFACC. Subroutine ACC also calculates THTPP(I), the second derivative of 
the function 

The local variables are defined as follows: 


ACCX 

EXACC 

KXPRI 

RELCOR 

SXACC 


Normalized acceleration of a disk 

Normalized acceleration of a disk due to rf electric field force and voltage- 
jump force 

Integer indicating location of table values for current value of | 

Relativistic correction factor 

Normalized acceleration of a disk due to space charge forces 
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XPRI Current value of | measured from beginning of current cavity 

XREL Current value of ^ measured from the last value of | for which space 

charge forces were calculated 


Subroutine APPLE 


Subroutine APPLE is called by the main program. Subroutine APPLE manipulates 
plot data and calls plot routines for plotting Applegate diagrams. 

The local variables are defined as follows; 


BREAK 

KEND 

NCURVS 

NPTS 

XDOWN 

YACROS 


Logical variable; true if two successive data points are more than 180° 
apart (indicates these two points should not be joined on plots by straight 
line segment) 

Logical variable; true if the current data point is the last data point 
Number of curves (NCURVS = 5 N^) 

Number of points per curve 
X array of points to be plotted 
Y array of points to be plotted 


Subroutine BTWNC 

A flow chart of subroutine BTWNC is given in figure 5. Subroutine BTWNC is 
called by subroutine BTWNS when the end of a cavity is reached. Subroutine BTWNC 
determines whether another pass throu^ the cavity is required. A second pass is not 
required if either of two conditions is met: (1) a second pass has already been made, or 
(2) the induced forward voltage is sufficiently close to the approximation of the induced 
forward voltage obtained by assuming constant disk velocities throughout the cavity. 

If a second pass is not required, the subroutine sets up for entrance into the new 
cavity. If needed, new tables are calculated by calling subroutine TABLE. Since a sec- 
ond pass through the new cavity may be needed, it is necessary to store the state at the 
beginning of the new cavity. Subroutine FSAPP is called to calculate DVTEMP, the ap- 
proximation to the induced. forward voltage for the new cavity. The new forward voltage 
is then calculated by phase shifting and attenuating the old forward voltage and adding 
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DVTEMP. If plots are desired, an approximation to the forward voltage for the cavity 
following the new cavity is required. We calculate this voltage in the same manner that 
we calculated the forward voltage for the new cavity. 

The local variables are defined as follows; 

KSWICH Determines which way to branch after calculating VGAP(KCAV) 

H or Appending H or HOLD to a variable name indicates that the variable is be- 

HOLD ing held in temporary storage in case it is needed again when a second 

pass through a cavity is made. 

T or Appending T or TEMP to a variable name indicates that the variable is be- 

TEMP ing held in temporary storage. 

Subroutine BTWNS 

Subroutine BTWNS is called by MAIN. Subroutine BTWNS determines what action 
should be taken between integration steps. First, the subroutine determines whether a 
new space charge calculation is required. A new calculation is required if the last cal- 
culation was KSCF steps ago or if the current value of | is at a boundary between two 
beam-diameter regions. Next, the subroutine determines whether the end of a cavity 
has been reached. If so, subroutine BTWNC is called. Finally, the subroutine deter- 
mines whether subroutine PLOTS should be called. Except in the first and last cavities, 
subroutine PLOTS is called four times per cavity. 

The local variables are defined as follows: 

LSCFIT Logical variable; true if current value of ^ is at a boundary between two 
beam- diameter regions 

NN(1) Beam diameter changes from first value to second value at end of cavity 

NN(1) 

NN(2) Beam diameter changes from second value to third value at end of cavity 

NN(2) 

Subroutine CAVP 

Subroutine CAVP is called by surboutine INDAT. Subroutine CAVP calculates 
quantities related to cavities and cavity strings. A cavity string is a string of con- 
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secutive cavities all having the same electrical and geometrical properties. The num- 
ber of cavities cannot exceed 650, and the number of cavity strings cannot exceed 35. 


Subroutine CKPDZ 

Subroutine CKPDZ is called by surroutine INDAT. Subroutine CKPDZ calculates 

th. 

KPDZC(K), the ratio of Pierce impedance to total impedance, for the k cavity. 


Subroutine CROOT 

Subroutine CROOT is called by subroutine SMSIG. Subroutine CROOT calculates 
the roots of the quartic polynomial involved in calculating small- signal parameters. 
The subroutine first solves for the roots of the cubic polynomial using the Newton- 
Raphson method. Using these roots as the first iterates, the subroutine then uses the 
Newton-Raphson method to solve for tiie roots of the quartic polynomial. 


Subroutine DVSUM 

Subroutine DVSUM is called by subroutine XINT. Subroutine DVSUM calculates one 
term in the summation for the induced forward voltage, DELV, and one term in the 
summation for the induced backward voltage, DELVR. These terms are added to the 
running sums for DELV and DELVR. The summations are those given in equations (33) 
and (34). 


Subroutine ENBAL 

Subroutine ENBAL is called by subroutine BTWNC. Subroutine ENBAL calculates 
the power balance at the end of the cavity in question. There are five terms in the 
power balance; PKE, PJUMP, PFW, PBW, and PLC. PKE is the kinetic power minus 
the initial kinetic power. PJUMP is the power due to all voltage jumps up to the cavity 
in question. PFW is the power in the forward wave minus the power in the forward 
wave in the first cavity. PBW is the power in the backward wave. PLC is the sum of 
power losses up to the present cavity. Power balance is equal to one plus the sum of 
these five terms divided by the initial beam power. 
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The local variables are defined as follows: 


ATTENF 

ATTENR 

ICL 

lOBVOB 

PBIP2 

PFl 

PFI 

PFIMl 


Attenuation factor for forward wave 
Attenuation factor for backward wave 
Number of the cavity in question 
Product of beam current and beam voltage 
Power in backward wave in cavity ICL + 2 
Power in forward wave in first cavity 
Power in forward wave in cavity ICL 
Power in forward wave in cavity ICI - 1 


Subroutine FSAPP 

Subroutine FSAPP is called by subroutine BTWNC. Subroutine FSAPP calculates 
an approximation, DVTEMP, to the induced forward voltage by assuming the disks have 
constant velocity in the cavity under consideration. 


Subroutine INDAT 

Subroutine INDAT is called by the main program. The subroutine first reads input 
data from a prestored data set. It then reads data entered at the terminal in NAMELIST 
format. Only input data that differ from those in the prestored data set need to be en- 
tered. The input data are stored in a data set for possible future use. The subroutine 
then calculates constants and initializes variables. The initial state is obtained in one 
of three ways: (1) if the input parameter ISTATE is 4, the initial state is read from a 
data set; (2) if ISTATE is not 4 and INITCV (initial cavity) is 1, the initial state is cal- 
culated from input data; (3) if ISTATE is not 4 and INITCV is not 1, the initial state is 
obtained from COMMON/STATE/. There is storage allotment for three states in 
common/state/. Which of the three is chosen is determined by ISTATE (see defini- 
tion of ISTATE in the section on description of input data) . 

After the initial state is obtained, it is stored in COMM ON/ST ATE/ in the first of 
the three storage allotments. This is done because the initial state of the present case 
may be used as the initial state in the following case, and thus it must be stored. The 
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initial state is also stored in the third of the three storage allotments if the input par- 
rameter ISAVE is an odd integer. This is done because the initial state of the present 
case may also be required as the initial state in some future case after the following 
case, and thus it must be stored. 

The input parameter ISAVE is the sum of II, 12, and 13 (see definition of ISAVE in 
the section on description of input data). If 12 is 2, the initial state of the present case 
is stored in a specified data set for possible use in a future case. 

If INITCV is 1, subroutine INDAT calculates the forward voltage for the first cav- 
ity. If INITCV is not 1, the forward voltage for the cavity previous to cavity INITCV is 
required for startup. This voltage will be known either from the previous case or by 
reading it from a data set. The backward voltages are calculated next. Finally, the 
induced backward voltages DELVR are stored in the array DVRTEM. This is done be- 
cause in the equations of motion we must use the induced backward voltages of the pre- 
vious case. We cannot use the DELVR array because this array changes during the 
simulation as the new induced backward voltages are calculated. We therefore use the 
DVRTEM array in the equations of motion. 


Subroutine NATTN 

Subroutine NATTN is called by subroutines' BTWNC and SCFIT. Subroutine NATTN 
calculates the number of the interception region in which XXX is located. 


Subroutine NR 

Subroutine NR is called by subroutines BTWNC and SCFIT. Subroutine NR calcu- 
lates the number of the beam-diameter region in which XXX is located. 


Subroutine OUTPT 

Subroutine OUTPT is called by subroutine BTWNC. Subroutine OUTPT calculates 
and prints output data. Since the output data have already been described in another 
section, we will not repeat the description here. 
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Subroutine PLOTS 


Subroutine PLOTS is called four times per cavity by subroutine BTWNS. For every 
other disk i, subroutine PLOTS calculates (p^^, the phase of the forward voltage seen 
by the i*'^ disk when it is at the middle of the gap in the k*"^ cavity. Althou^ there is 
only one per cavity per disk, it is desirable to plot more points per cavity per disk 
to obtain smoother plots. An interpolation procedure is used to obtain four points per 
cavity per disk. When | is in the first half of the cavity, the subroutine calculates a 
voltage W by interpolating between voltages in the previous cavity and the present cav- 
ity. When I is in the second half of the cavity, the subroutine calculates W by inter- 
polating between voltages in the present and following cavities. The phase angle is 
then calculated as the phase of VV. 

The local variables are defined as follows: 

PHIK Phase angle <p^^ 

SHALF Logical variable; true if | is located in second half of current cavity or if 
current cavity is cavity INITCV 

T1 Normalized time at which disk is at center of previous gap 

T2 Normalized time at which disk is at center of present gap 

T3 Normalized time at which disk is at center of following gap 

VI Forward voltage, in previous cavity, whose phase is that seen by disk when 

it is at center of gap 

V2 Forward voltage, in present cavity, whose phase is that seen by disk when 

it is at center of gap 

V3 Forward voltage, in following cavity, whose phase is that seen by disk when 

it is at center of gap 

VV Interpolated value of forward voltage 

XI Normalized axial position at center of previous gap 

X2 Normalized axial position at center of present gap 

X3 Normalized axial position at center of following gap 
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Subroutine RFACC 


Subroutine RFACC is called by subroutine ACC. Subroutine RFACC calculates 
EXACC, the normalized acceleration of a disk due to rf electric field force and voltage- 
jump force. 

The local variables are defined as follows: 

EXACC Normalized acceleration of a disk due to rf electric field force and voltage- 
jump force 

VF Forward voltage at time disk arrives at current value of ^ 

VR Backward voltage at time disk arrives at current value of ^ 

*• " • ■ ■■ ‘ 

*'.♦ Subroutine SCACC r' 

Subroutine SCACC is called by subroutine SCFIT. Subroutine SCACC calculates the 
space charge acceleration of one disk due to another by linear interpolation on space 
charge force tables. Although we use the term space charge force tables, these tables 
actually evaluate accelerations rather than forces. If the two disks are in two different 
beam- diameter regions, two linear interpolations are done and the average is taken. 

The local variables are defined as follows: 

NREF Number of beam-diameter region of reference disk 

NSRC Number of beam-diameter region of source disk 

SCAO Space charge acceleration of reference disk ; 


Subroutine SCAT 

Subroutine SCAT is called by subroutine INDAT. For each beam diameter IB and 
th 

for the I ” separation distance X, subroutine SCAT calculates the space charge accel- 
eration SCATMP (I, IB) of one disk due to another. Using EQUIVALENCE statements, 
we put SCAl, SCA2, and SC A3 into COMMON/MREAL/ instead of SCATMP. The 
arrays SCAl, SCA2, and SCA3 make up the space charge force tables. In this sub- 
routine, X is used as the separation distance between disks. Since X is also used as a 
global variable (current value of |), it is necessary to save the global variable in 
XSAVE. At the end of the subroutine, the value of the global variable is returned to X. 
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Subroutine SCAT also calculates ZRATIO, the ratio of the largest space charge acceler- 
ation to the smallest. The smallest space charge acceleration occurs when the separa- 
tion distance is largest. The largest separation distance is determined by the input 
variable SCMAX. 

The local variables are defined as follows: 


LDISKR 

SCAMAX 

SCATMP 
a, IB) 

X 

XSAVE 

ZMAX 

ZRATIO 


Relativistic thickness of a disk, m 

Maximum space charge acceleration of one disk due to another disk 

For the IB*"^ beam diameter and the separation distance, the space 
_ charge acceleration of one disk due to another 

Separation distance between disks 

Temporary storage for global variable X 

For the largest separation distance, ZMAX is the largest space charge 
acceleration among the beam diameters 

The ratio of the largest space charge acceleration to the smallest 


Subroutine SCFIT 

Subroutine SCFIT is called by subroutine BTWNS. Subroutine SCFIT first calcu- 
lates what the space charge accelerations would be if the present curve fits are used. 
Then for each disk the subroutine calculates the space charge acceleration of the disk 
by summing the accelerations due to all the other disks. The acceleration of one disk 
due to another is obtained by linear interpolation on the space charge force tables. The 
coefficients of the new quadratic curve fits are then calculated. Finally, the subroutine 
determines whether the frequency of the space charge calculation should be changed. 

This is done by comparing the space charge accelerations based on the old curve fits 
with those just calculated. If the difference for any disk is greater than the product of 
TOLSC and the maximum acceleration due to the rf electric field in the cavity, the fre- 
quency of calculating space charge accelerations is doubled. 

The local variables are defined as follows: 

ARG Normalized separation distance between reference disk and source disk 

ERRSC(I) For the i^^ disk, the difference between the space charge acceleration based 

on the old curve fits and the space charge accelerations just calculated 
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NBSRC 

NREF 

SCAMAX 

SXACC(I) 

THTREF 

XREL 


Number of beam-diameter region of source disk 
Number of beam- diameter region of reference disk 
Maximum value in array ERRSC 

For i^'” disk, the space charge acceleration based on the old curve fits 
Current value of 9 for reference disk 

Current value of | relative to last value of ^ for which space charge 
forces were calculated 


Subroutine SMSIG 

Subroutine SMSIG is called by the main program. Subroutine SMSIG calculates 
small- signal parameters as given in reference 8. The Pierce C parameter is given by 



where K 

P 


is the Pierce impedance. 


The Pierce velocity parameter is given by 



coC 




a^' 

The Pierce loss parameter is given by 


d = 


oL 

20 


logfoe 


cuLC 

^0 


The plasma frequency 




is given by 
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r>2xT 2 

R Ne 


“q V m_e 


e 0 


where R is the plasma frequency reduction factor obtained by table lookup and N is 


N =. 


e 


The Pierce Q parameter is given by 


Q = 




4C^(co 


The roots Sg, and 5g are those of the cubic equation 

1 


62 =. 


- 4QC 


-b + jd + j6 

The roots of the cubic equation are only approximate. The exact roots are those of the 
quartic equation 


g2 _ (l+Cb-jCd)(l+Cbr _ 4QC 
r<r •> ol 

-b + jd 


jd + jd - — 6^ + (b - jd)2 
2 - — 


The launching loss is given by 


+ Ag = 20 log 


6" + 4QC 


(5i - dgUdj^ - 6g) 


The gain per cavity in dB is given by 


54. 6 /^e (6JCL 

AG = i 

A 
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where \ is the beam wavelength and L is the length of the cavity. 

The local variables are defined as follows: 

Launching loss, A1 + A2 

Pierce velocity parameter b 

Pierce parameter C 

Pierce loss parameter d 

root of quartic polynomial 

Gain per cavity, dB per cavity 

Pierce impedance, ohms 

Pierce parameter Q 

Plasma frequency reduction factor 

SSGAIN(K) Small- signal gain at cavity, dB 

WQ Plasma frequency, rad/ sec 

Subroutine STDAT 

Subroutine STDAT is called by the main program. Subroutine STDAT calculates 
standard data, that is, constants whose values are independent of input data. 


A1PA2 

BB 

CC 

DDD 

DEL (I) 

DGAIN 

KP 

QQ 

R 


Subroutine TABLE 

A flow chart for subroutine TABLE is given in figure 6. Subroutine TABLE is 
called by subroutine BTWNC when a new cavity is entered and tables have to be com- 
puted for the new cavity. If the new cavity has the same properties as the old cavity, 
subroutine TABLE is not called. The subroutine calculates the following tables; 

TXR(n, 1), Real and imaginary parts of the table Snk needed to calculate induced volt- 
TXI(n, 1) ages; n = 1, . . . , NXGRID 
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TXR(n,2), Real and imaginary parts of the table needed to calculate rf electric 

TXI(n,2) field forces; n = 1, . . . , NXGRID 

TXR(n, 3) Table R^ needed to calculate voltage-jump forces; n = 1, NXGRID 
Prom equation (32), dropping the subscript k, we obtain 


TXR(n, 1) + i TXI(n, 1) = S* 



(39) 


From equation (21), dropping the subscript k, we can define such that 


TXR(n, 2) + i TXI(n, 2 ) = = 



(40) 


Prom equations (23) to (25), dropping the subscript k, we can define U^, and 

Y__ such that 
m 


TXR(n, 3) = R^ = Tq + cos^^^^ (41) 

m=l 


or 


TXR(n, 3) = R^ = Uq + 


00 



m=l 


(42) 


Equations (39) and (40) can be rewritten as 
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TXR(n, 2) + i TXI(n, 2) = e 


-%yn 



m=- oo 


(43) 


(44) 


From equations (41) to (44) , it is clear that a significant amount of the computer time 
for calculating tables is spent on evaluating cos(27rmyj^/L) and sin(27rmy^/L) . Sub- 
routine TABLE uses a scheme that significantly reduces the computer time for eval- 
uating these trigonometric functions. For given values of m and n, we define 
and k2 such that 



(2n - 1) = k^(2N^) + kg 


0 < kg £ 2N^, - 1 


Then it can be shown that 


cos 




0 < kg < 2N^ - 1 


sin 




0 ^ kg ^ 2N^ - 1 


Thus, the trigonometric functions take on only 2N distinct values. These values are 
calculated once in subroutine INDAT and stored in the arrays COSY and SINY. Sub- 
routine TABLE needs only to calculate kg to obtain the trigonometric functions. 

Subroutine TABLE calculates the tables for positive values of y^ (ri^t half of cav- 
ity) . Symmetry is then used to obtain the tables for negative values of y^ (left half of 
cavity) . 
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The local variables are defined as follows: 


COSBOY(N) 
COSDEL 
FAC(1,NN) 
FAC(2,NN) 
FAC(3,NN) 
FAC (4, NN) 
FAC(5, NN) 
LJUMP 


cos()3,qZ^), where k is the number of the cavity in question 
cosi^oL^'i^^), where k is the number of the cavity in question 
Array used in calculating TXR(IZ, 1) and TXI(1Z, 1) 

Array used in calculating TXR(IZ, 2) and TXI(IZ, 2) 

Array used in calculating TXR(IZ, 3) 

Array used in calculating TXR(IZ, 3) 

Array used in calculating TXR(IZ, 3) 

Equals zero if TXR(IZ, 3) is to be calculated; equals 1 otherwise 


NMAXPl 

NNMAX 

NXDUBL 

NXHALF 

SINBOY(N) 

SINDEL 

Y 


NMAX+1 

2*NMAX+1 

2*NXGRID 

NXGRID/2 

sin(^QZ^), where k is the number of the cavity in question 
sin(/3QLj^/N^), where k is the number of the cavity in question 
Axial position relative to center of gap 


Subroutine XINT 

Subroutine XINT is called by the main program. Subroutine XINT integrates the 
equations of motion over one step. The numerical integration method is described in 
the section on the dynamics of beam disks. 


CONCLUDING REMARKS 

The use of the coupled- cavity traveling wave tube for space communications had led 
to an increased interest in improving the efficiency of the basic interaction process in 
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these devices throu^ velocity re synchronization and ottier methods. In order to ana- 
lyze these methods, we have recently developed a flexible large- signal computer pro- 
gram for use on the IBM 360/67 time-sharing system. The present report is a users’ 
manual for this program. The report describes the program in sufficient detail to allow 
a user to make modifications in the program if desired. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 5, 1977, 

506-20. 
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APPENDIX A 


SYMBOLS 


a 

b 

^mk 

c 


e 



f 


^0 

J^(z) 

k 

m 

e 

Nb 

Nc 

Nd 

Nz 

Pk 

‘^nk 




^nk 

®nk 


t^(z) 

“O 


tunnel radius 
beam radius 

factor used in tables, defined in eq. (17) 

speed of light 

charge of electron 

axial force on i*"^ disk 

frequency 

beam current at tube entrance 

fundamental Fourier component of current density at 
til 

length of k ” cavity 

axial thickness of a disk 

half-length of k*^ gap 

mass of a disk 

mass of an electron 

number of beam- diameter regions 

number of cavities 

number of disks per beam wavelength 
number of nodes per cavity 
power, k cavity 
table for calculation of rf forces 
charge of a disk 

table for calculation of voltage-jump forces 
table for calculation of induced voltages 
reciprocal of frequency 

time of arrival of i^^ disk at axial location z 
initial beam velocity 



V, 


bk 


V, 


fk 


AVi 


bk 


AV. 


fk 


Vo 

V 


Jk 

Vj(z) 

^nk 

^nk 

AZk 

(“L)bk 

(^L)fk 

^mk 

^k 


OS^L), 


y 


mk 


n 


'■0 

0i(l) 


th. 

backward voltage, k*^” cavity 
th 

forward voltage, k''^ cavity 


th 

backward induced voltage, k*"" cavity 
forward induced voltage, k*"^ cavity 
beam voltage 

th 

voltage jump, k*^ cavity 
th 

gap voltage, k*' cavity 

fVi 

velocity of i''" disk at axial location z 

axial position, relative to center of k^^ cavity, of n^^ node in cavity 

th 

interaction impedance for k*"” cavity 

axial position, relative to beginning of first cavity, of n^^ node in k^^ 
integration step size in cavity 

th 

loss factor for backward voltage in k*^^ cavity 

th 

loss factor for forward voltage in k^“ cavity 
propagation factor, defined in eq. (18) 

fH 

phase shift of voltage for k cavity 
factor used in tables, defined in eq. (19) 
n*'^ coefficient for Simpson’s rule integration 
permittivity of free space 


•th 


disk at normalized axial location ^ 


normalized time of arrival of i 
wavelength 

p shaping parameter for electric field 

I normalized axial position; independent variable for equations of motion 

ct) angular frequency 

Subscripts: 

k cavity number 

mid middle 


cavity 
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Figure 3. - Flow chart for computational procedure. 


















Figure 4. - Applegate diagram. 
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Figure 5. - Flow chart for subroutine BTWNC. 
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Figure 6. - Concluded. 
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APPENDIX B 


SAMPLE PROGRAM RUN 

A sample run, which illustrates the use of the program, is given in this appendix. 
The output of this run is given in figure 7. The input parameters are chosen to repre- 
sent a tube with a lumped sever, a single-step circuit velocity reduction, and an accel- 
erating voltage jump. This run also illustrates the use of the iterative procedure to 
calculate the effect of the backward wave on tube performance. This run is not intended 
to represent a finished tube design but to illustrate use of the various program capa- 
bilities. 

A list of input parameters for the first case is followed by the output for the first 
case - a first iteration through the entire 37-cavity tube, ignoring the backward wave. 
Note that the tube is severed at cavities 13 and 14. Going from cavity 25 to cavity 26, 
we encounter a step reduction in circuit velocity. A voltage jump occurs at cavity 32. 
Figure 8 shows the Applegate diagram for this case. 

The remaining nine cases illustrate the iterative procedure used to evaluate the 
effect of the backward wave. The output of each case is preceded by a list of those in- 
put parameters that are changed from the previous case. 

Effectively, there are five iterations through the input section and five through the 
output section. Note that satisfactory convergence is achieved with three iterations in 
each section. Comparing the first case to the last shows the effect of the backward 
wave. A list of normalized final velocities of the disks in a beam wavelength is also 
included in the output for each case. 
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Figure 7. - Continued. 
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NORMALIZED VELOCITIES 

1.11805 1.08884 1.06199 1.03523 1.05514 1.05805 1.03391 0.98203 0.99133 0.99642 0.98320 0.97958 
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NOBHSLIZED VELOCITIES 

1. 11769 1.08871 1.06170 1.03506 1.05526 1.05811 1.03324 0.98214 0.99126 0.99719 0,. 98395 0.97967 

0.98538 0.99306 0.94213 0.85812 0.75699 0.71273 0.72882 0.74377 0^ 71353 0.71977 Oi 976«8 1.051.91 
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